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Abstract 



Finite-size scaling expressions for the current near the continuous phase tran- 
sition, and for the local density near the first-order transition, are found in the 
steady state of the one-dimensional fully asymmetric simple-exclusion process 
(FASEP) with open boundaries and discrete-time dynamics. The correspond- 
ing finite-size scaling variables are identified as the ratio of the chain length 
to the localization length of the relevant domain wall. 
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I. INTRODUCTION 



We consider the fully asymmetric simple-exclusion process (FASEP) on a finite chain of 
L sites with open boundaries. For mathematical definition of the exclusion processes we 
refer the reader to the book [|l|, and for a recent review on the relevant class of exactly 
solvable models for many-body systems far from equilibrium to [Q, see also 0]. We recall 
that each site i G {1,2, ...,L} of the chain is either empty or occupied by exactly one 
particle. The particles obey a discrete-time stochastic dynamics according to which they 
hop with probability p only to empty nearest-neighbor sites to the right. The open boundary 
conditions imply that at each time step (update of the whole chain) a particle is injected 
with probability a at the left end of the chain {i = 1), and removed with probability (3 at 
the right end {i = L). The order in which the local hopping, injection and particle removal 
takes place is specified by one of the basic discrete-time updates, see Here we explicitly 
consider the case of forward-ordered sequential update. 

We mention that the case of random-sequential update was solved first by using the re- 
cursion relations method 0,0, and then by means of the matrix-product Ansatz (MPA) [0 . 
Next, the method of MPA was successfully applied for obtaining the steady-state properties 
in all the basic cases of true discrete-time dynamics: forward- (— and backward-ordered 
(^) sequential 0, [|T^], sublattice-parallel (s- ||) [|TT|, |T^, and, finally, fully parallel 
dynamics [jl3[ , [1^ . 



The phase diagram for all the discrete-time updates has the same structure, as shown 
in Fig. 1: it contains a maximum-current (m.c), low-density (l.d.), and high-density (h.d.) 
phases. The maximum-current phase is separated by lines of continuous phase transitions, 
(a = ttc, /9c < < 1) and (/5 = Pc, < a < 1), from the low-density and high-density 
phases, respectively. Here etc and Pc are the critical values of the injection and removal 
probabilities: 



a, = (3c = l-^Jl~p. (1) 

The above phases were identified with respect to the analytic form of the bulk current: for 
fixed p, the current in the low- (high-) density phase depends only on a (/3), and in the 
maximum-current phase it is independent of both a and (3. On crossing the borderline 
between the maximum-current and the low- (high-) density phase, the current itself and 
its first derivative with respect to a [jS) change continously, and the second derivative with 
respect to a {/3) udergoes a finite jump. The coexistence line between the low- and high- 
density phases is given by (a = /5, < (3 < /3c); on crossing it the bulk density undergoes a 
finite jump. 

Here, the exact finite-size expressions for the current and local density in the steady 
state of FASEP with open boundaries and forward-ordered sequential update, derived in 
IIIOII , are analysed within the framework of finite-size scaling (FSS) at continuous (for the 
current) and first-order (for the density) phase transitions. The appropriate scaling variables 
are identified and the corresponding scaling functions for the current and local desity are 
explicitly obtained. The notions of the Privman-Fisher anisotropic FSS have been recently 
extended to non-equilibrium systems belonging to the directed percolation and diffusion- 
annihilation universality classes, f^, ||15|, [|16|. To the best of our knowledge, the present 



study is the first step in the analytic confirmation of FSS for an exactly solved model of 
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a driven lattice gas with open boundaries. Since we are dealing with phase transitions in 
the steady-state of the FASEP with discrete-time updates, the equvalent two-dimensional 
lattice model is infinite in the temporal direction but finite in the spatial one. 

Note that the current has the same value in the cases of forward-ordered (^), backward- 
ordered (<— ) sequential, and sublattice-parallel (s- ||) dynamics, i.e., = J|7 = ". The 
corresponding local densities at site i G {1, . . . , L} are related to each other 

As shown in [|l3l, the current j| and local density for the FASEP with fully parallel 
update can be simply expressed in terms of those for the forward-ordered sequential update: 

4^^. plw^^fliS. (3) 

Due to these relations, the results derived here suffice to explicitly obtain the current and 
density FSS functions for each of the basic discrete-time updates. 

Concerning the notation, we note that the exact finite-chain results obtained in [|10| are 
conveniently expressed in terms of the parametrs 

p — a P — P 



d=sjl-p, a = d + d-\ e = ^, V = ^, (4) 



which will be used here too. 



II. FINITE-SIZE SCALING AT THE CONTINUOUS PHASE TRANSITION 



Let us consider first the continuous phase transition across the boundary a = Oc, (3c < 
(3 < 1 between the low-density phase and the maximum-current phase. In terms of the 
variables (^) the equation of this boundary reads ^ = 1, —d < r] < 1; the m.c. phase 
occupies the region —d < ^ < 1, —d < r] < 1, and the region ^ > I > r] > —d, called region 
All, lies in the l.d. phase, see Fig. 1. 

According to the basic FSS hypotheses, the FSS variable in the case of a continuous 
transition, characterized by diverging bulk correlation length A, should be given by the ratio 
L/X, where L is the finite-size of the system. As it is well known, in the case at hand the 
inverse correlation length A^^ in the l.d. phase is, see, e.g., |T^, |T^, 



A^^ = In 



1 + 



e(a + 2) 



(e>l) 



and A^ ^ = in the m.c. phase. Hence, the FSS variable as L — > oo and ^ 
to be given by 



(5) 

1+ is expected 
(6) 



where 
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Xi 



Ci = {a + 2)-^/\ t = ^-l. 



(7) 



We emphasize that here we study a boundary induced non-equihbrium phase transition, 
and the physical quantity which measures the distance from the steady-state critical point 
is related to the injection probability: t = ^ — 1~1 — 

Consider now the finite-size current JL{^,ri) = Zi_i{^,vi) / Zl{^,vi), where Zi^^^r]) is the 
normalization constant of the steady-state probability for a lattice of L sites. The following 
exact representation of Zi{^,ri) in the subregion All of the low-density phase > 1 > ?7) 
has been found in |[T0[| : 



.P. 



(8) 



Here the expression for the normalization constant in the maximum- current phase 7^ 

L r 



1 



(9) 



involves the integral 



vr JO 1 



(a + 2 cos^)"^ sin^ 



- 2^ COS(/) + ^2 



(10) 



which is a non-analytic function of at ^ = 1. For all finite L the normalization constant 
Z^^^(^,?7) in region All represents an analytic continuation of Zf^'^'{$^,r]) from the domain 
1^1 < 1 to the domain |^| > 1, see jlOl . 

A direct application of Laplace's method for evaluation of the integral (|1^) as L ^ oo 
shows that it changes its leading-order asymptotic behavior from 0{L~^^'^) for ^ 7^ 1 to 
0{L-^/'^) for ^ = 1, see Eq. (14) in [0. The finite -size expression that interpolates between 
these limiting asymptotic forms can be readily obtained by using small- argument expansion 
of the trigonometric functions in the integrand. The result for Xi{L,t) = 0(1) is 



(a + 2) 



L+1/2 



-X(\x,\ 



(11) 



where the FSS function X(-) is given by 

X{x) = l- v^xe^'n - $(x)], $(x) = ^ 

vvr Jo 



dt. 



Thus, keeping the 1/L corrections to the finite-size scaling form, we obtain 



(12) 



zri^,v) 



.p. 



(a + 2) 



L+1/2 



X(|xi| 



7] a + 2 
(l-r/)2 2L 



(13) 



The small- and large-argument asymptotic behavior of X(-) readily follows from that of the 
Fresnel integral $(■): 
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X{x) 



1- ^x + Oix^), 

- + 0(x-6) 



as X — > 0, 
as X ^ oo. 



(14) 



Let us first evaluate the asymptotic behvior of tlie finite-size correction to tlie bulk 
current in the maximum-current phase, 



m.c. 

oo 



p 



(15) 



d{a + 2y 

as L oo and ^ — 1~ at fixed t] < 1. We readily obtain in the leading order of magnitude, 



^[(a + 2)h-iiO-lLm 



|l-e|3/2 2 



dx- 



^3/2 vr Jo 1 + x^ 
[l-2x?X(|xi|; 



(a + 2)3/2 



Hence we derive 



2v/?L3/2 



p 



L 2d{a + 2) 



2x2X(|xi|) 



X(|xi|) 



(16) 



(17) 



The asymptotic behavior of the above finite-size correction to the current as Xi — (say, 
as ^ — *• 1" at large fixed L) follows from Eq. ([T^): 



p 



L2d{a + 2)' 



(18) 



In the limiting case xi — — oo (say, L — > oo at ^ < 1 fixed close to unity), we obtain that 
the finite-size correction is again of the order 0{L~^): 



3p 



L2d{a + 2)' 



(19) 



Consider next the asymptotic behvior of the finite-size correction to the bulk current in 
the low-density phase. 



^r(e,r])-jr-(o 



p 



as L — ^ oo and ^ — > 1+ at fixed r] < 1. With the aid of the identities: 

e[(a + i + r')/L~i(o - hm = h-i{o), 

we obtain in the leading order: 



l.d. 



P 



L2d{a + 2) 



1 



20Fxie^i +X(xi)_ 



(20) 



(21) 



(22) 



(23) 
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Hence, by using Eq. (jTj) we readily derive the asymptotic behavior of the above finite-size 
correction — s> 0^ (say, as ^ ^ 1^ at large fixed L): 

(24) 

In the other limiting case xi ^ oo (say, L — oo at ^ > 1 fixed close to unity), we obtain 
that the finite-size corrections to the current are exponentially small: 

Jt'\^,v)-Jl'iO = 0(x^'e-^'A. (25) 



As a finite-size order parameter of the continuous phase transition one could consider 
the difference in the finite-size currents: 

AL{^,r^):=jri^,v)-Jt'\^,v) {^>l,V<l)- (26) 

Here </™''^ (^, t]) represents the analytic expression for the finite-size current in the maximum- 
current phase, i.e., J™'^' = Zf^_f{/Zf^-'^- with Zf^-^- defined by Eqs. (^ and (0), evaluated 
in subregion AH of the low-density phase (at ^ > 1, ?7 < 1). For the corresponding bulk 
quantity we have 

Aoo(0 = - ^'-(0 = %f^, (27) 

p{l - a) 

which suggests the critical exponent for the order parameter (3 = 2. In the finite-size case, 
by taking into account that 

P P ^ 1 P 4 /28^ 

d{a + 2) d(a + ^ + ^-i) L d{a + 2) 1 + xj/ L' ^ ' 

and combining the above results, we obtain in the leading order 

^"^^'''^ ^ L2d(J+2) {x(xO [xi + 4^e-?X(xO] } ' ^^^^ 

Since the high-density phase maps onto the low-density phase under the exchange of 
arguments ^ rj (equivalently, a ^ [3), the FSS properties of the continuous transition 
across the boundary (3 = (3c, between the high-density phase and the maximum- 

current phase follow trivially from the above results and the particle-hole symmetry. 



III. FINITE-SIZE SCALING AT THE FIRST-ORDER TRANSITION 

In the thermodynamic limit the first-order phase transition, which occurs across the 
borderline /? = a, < a < (^7 = ^ ^ 1) between subregions AI and BI, manifests itself 
by a finite jump in the bulk density: 

<;1.(^)-p1;^i.(0|^,^ = J;f^^, >o- (30) 
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Quite peculiarly, this transition is characterised by another diverging correlation length 



^ ' a + 7] + 7] 
which suggests the finite-size scaling variable 



-\^-v\ + oi\^-v\ 



£2 := L/X ~ C2\h\L, C2 



1—7] 



-2 



h = ^-7]. 



(31) 



(32) 



a + 7] + 7]^^ 

Our analysis starts with the exact result found for the finite-chain local density p£(i; ^, 7]) 



in region D=AIUBI of the phase diagram, see Eqs. (4.21) and (A13) in [10|. This result can 
be cast in the form (^ 7^ ?]): 



p?(^;e,^) 



p 



(33) 



where. 



'd\ (a + 2) 



L-l 



,PJ (^ - v)Zl L 

_ fd\'' (a + 2)^-1 



/ri\L-i vm.c. 



(a + 2) 



L-i 



/■n\»-l 7m. c. 



+ 



d 1 



d/ (a + 2)^-1 



(34) 



For brevity of notation we have omitted the arguments ^ and rj of the function Zl{^,7]). 
Here the term FL{i) = FL{i;^,7]) is the antisymmetric (with respect to the center of the 
chain) function of the integer coordinate i, FL{i;^,7]) = —Fl{L — i + 1;^,?]), defined for 
1 < i < [L/2] (where [x] denotes the integer part of x) by the equation 



L-l 



\P/ n=0 

The normalization constant in the region D=AIUBI > 1,7] > 1) for ^ 7] is given by 



(35) 



z^{^,v) 



'dy (a + 2) 



(e - r')e^/'« -iv- r^)e^/^^l + Zn^, T]). 



yPj 

Let us analyse this expression when = 7] + h, h ^ 0. For Z'^'^'{7] + h, 7]) we have 



(36) 



zriv + h,v) 



[{7] + h)h{7] + h)- 7]Il{7])] 



[p) h 



(37) 



where 



'a + 2 cos0)^ sin^ 



TT J (1 — 2?7 COS + r^^)^ 
For ^ = T] + h, asL^oo and h ^ we have 

L/\+h = L/X, + C2hL + 0{Lh^), 



hence 



{a + 2Y 



pC2hL _ 1 



(38) 



(39) 



(40) 



Next, by using the upper bound hiO < (a + 2)^/o(0) where /o(0 = 1 ^^i \^\ < 1 and 
loiO = 1^1 — 1) easily obtains that in region D 



\Fl{i;^,v)\ <-- 



.P. 



L-l 



(L- l)(a + 2) 



L~2 



(41) 



Therefore, in view of Eq. pO|), the contribution of the term proportional to Fi{i;^,ri) into 
the local density is exponentially small. Thus, only the first three terms in the right-hand 
side of Eq. ( P^ contribute into the leading-order expression for Cl^{ri + h,ri): 



a + rj + rj 



„C2hi _ 1 
- + (^-^ \c2hL_l 



(42) 



In deriving the above expression we have assumed that both z ^ 1 and L — i ^ 1. Finally, 
by taking into account that 



J^{n + h,r^) 



P 



d{a + T] + 7] 



—,+0{h), 



(43) 



we obtain the leading order expression for the local density on the macroscopic scale i/L = r, 
< r < 1: 



1 



a + 7] + 7] ^ 



d + T] ^ + {t] - 7] 



^X2r 



qX2 _ I 



(44) 



Here we have introduced the FSS variable X2 = C2hL, compare with Eq. (p2D . In the 
limit X2 —>■ the above expression reduces to the well-known linear density profile on the 
coexistence line ^ = rj: 



pLirL;V,v) 



d + T] ^ -\- [r] — 7] 



a + r] + 1] 



-1 



(45) 



In the limit X2 —>■ +00 (x2 —>■ —00) one recovers the bulk density in the low-density (high- 
density) phase. 



8 



IV. DISCUSSION 



The mathematical mechanism of the phase transitions in the FASEP with open bound- 



aries has been revealed in [10| as qualitative changes in the spectrum of the lattice translation 



operator C ||T8[. In the region ^ < 1, ?7 < 1, occupied by the maximum-current phase, the 
spectrum is continuous and fills with uniform density the interval from {d/p){a — 2) to 
{d/p){a + 2). When f] < 1 but ^ becomes larger than unity, an eigenvalue 



(46) 



splits up from the continuous spectrum and dominates the properties of the low-density 
phase in region AIL Due to particle-hole symmetry, when ^ < 1 but t] exceeds unity, the 
eigenvalue that splits up from the continuous spectrum and dominates the properties of the 
high-density phase in region BII is z/(?7). Thus, the logarithm of the ratio of the corresponding 
eigenvalue to the ceiling of the continuous spectrum defines the relevant inverse correlation 
length A^'^ or A~^, see Eq. (|). When both ^ > 1 and r/ > 1, which is the case in 
region D=AIUBI, there are two eigenvalues, z/(^) and i^(?7), above the continuous spectrum. 
Obviously, these eigenvalues become degenerate on the coexistence line ^ = rj > 1, which 
explains the appearance of the diverging correlation length (^). Note that the explicit 
expressions for the correlation lengths depends on the type of update only: they are the 
same for all the true discrete-time updates, and for the random sequential update, see Eq. 
(78) in [0, 



A 



-1 



In 



;i/2 - a) 



a(l — a) 



(47) 



The physical meaning of the above correlation lengths emerges in the domain wall picture 
developed in [r^|: the lengths Ag, A^, and A are interpreted as localization lengths of the 
domain walls between the low-density/maximum-current, high-density/maximum-current, 
and low-density/high-density phases, respectively. The complete delocalization of the low- 
density /high- density domain wall on the coexistence line explains the linear density profile 
(^): it is the result of ensemble averaging over configurations with uniform probability 
distribution of the domain wall position 0. 

Thus, the FSS variable for any of the phase transitions in the FASEP with open bound- 
aries has the physical meaning of a ratio of the chain length to the localization length of the 
relevant domain wall. The explicit expressions for the FSS functions have been derived here 
for the model with forward-ordered sequential update. The corresponding expressions for the 
other basic discrete-time updates follow under the mappings mentioned in the Introduction. 
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FIG. 1. The phase diagram in the plane of the injection and removal probabilities a and /? (see 
the text) for hopping probability p=0.75. The maximum-current phase occupies region C. Region 
A=AIUAII corresponds to the low-density phase, and region B=BIUBII to the high-density phase. 
Subregions Al (Bl) and All (BII) arc distinguished by the different analytic form of the density 
profile. The boundary between them, /3 = /3c, < a < ac (a = ac, < /9 < /3c), is shown by 
dashed segment of a straight line. The solid line a = (3 between subregions Al and Bl is the 
coexistence line of the low- and high-density phases. The curved dashed line is the mean-field line 
(l-a)(l -/?) = ! -p. 



